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MAP  OF  AN  INVERSE  PROBLEM  FOR  THE  WAVE 
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Department  of  Mathematical  Sciences 
Rice  University 
Houston,  Texas  77251-1892 


1  Introduction 

A  simplified  model  which  governs  many  physical  processes  such  as  seismic  and  acoustic  wave 
propagation  is  the  following  linear  acoustic  wave  equation: 

=  (1.1) 

where  a  —  is  the  logarithm  of  the  density,  c  =  c(x )  is  the  sound  speed  of  the  medium, 
and  /  =  f(x,t)  is  the  source  term  which  introduces  the  energy  to  the  problem.  If  er,  c 
and  /  are  given  along  with  appropriate  side  conditions,  the  forward  (or  direct)  problem  is 
to  determine  u  =  u(x,t.),  the  excess  pressure.  For  appropriate  choices  of  a,  c,  and  /,  u  is 
determined  uniquely  by  standard  linear  hyperbolic  theory  of  partial  differential  equations 
(p.cl.e.).  Thus  the  problem  stated  above  defines  a  map  from  the  coefficients  to  the  solution 
of  the  wave  equation.  In  this  paper,  we  study  an  aspect  of  the  regularity  of  this  map,  or 
rather  its  composition  with  the  trace  on  a  time-like  hypersurface. 

Throughout  this  work  we  shall  restrict  ourselves  to  the  special  case  of  constant  velocity 
c,  though  we  believe  that  the  ideas  in  this  work  may  be  extended  to  cover  some  more  general 
cases. 

To  fix  ideas,  write  x  E  IRn  as  (x\  ,r„),  where  x'  E  IR"-1 ,  xn  €  IR.  We  assume  that  the 
problem  is  set  in  thfe  whole  space  IRn  and  u  =  0  in  the  past  (t  <  0).  Take  f(x,t)  =  6(x,t) 
as  an  ideal  point  source.  Thus  u  is  the  retarded  fundamental  solution: 

n«  —  V<7  ■  Vi/.  =  6(x,t) ,  (x,t)  E  Htn  x  IR  (1-2) 

M  =  0  ,  t  <  0  ,  (1.3) 

where  □  is  defined  to  be  d\  —  A,  and  A  is  the  Laplacian. 

1  This  work  was  partially  supported  by  the  National  Science  Foundation  under  grant  DMS  86-03614  and 
DMS  89-05878,  by  the  Office  of  Naval  Research  under  contracts  N00014-K-85-0725  and  N00014-J-89-1 115, 
by  AFOSR  89-0363,  and  by  the  Geophysical  Parallel  Computation  Project  (State  of  Texas). 
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Define  the  forward  map  F  as: 


F  :  a  -  («  |In=„  ,  (1.4) 

where  f>  G  C'£°(IRn+1)  is  supported  inside  the  conoid  {f  >  |x|}  and  near{a:„  =  0}. 

Because  F  is  nonlinear,  one  wants  to  work  with  the  formal  linearization  (or  formal  deriva¬ 
tive)  DF ,  with  respect  to  the  reference  state  (<to,  u0),  defined  by  first  order  perturbation 
theory  (Born- approximation ) . 

OSu  —  Vao  ■  VSu  =  VSa  •  Vu0 
8u  =  0  ,  t  <  0  . 

The  formal  derivative  DF(ao)  is  given  bj? 

DF(a0)6a  =  (<j>8u)  \Xn=0  .  (1.7) 

It  is  our  main  goal  in  this  work  to  determine  appropriate  spaces  of  the  domain  and  range  of 
F  for  which 

the  formal  derivative  DF  is  bounded. 

Throughout  this  paper,  we  shall  always  assume  that 

(A)  supp(8c t)  C  {xn  >  e}  . 

for  e  >  0  small.  In  some  applications,  this  assumption  is  realistic,  as  the  density  can 
be  measured  directly,  near  the  location  of  receivers  (i.e.  xn  —  0).  Also  the  domain  of 
dependence  properties  of  the  wave  equation  implies  that  <f>u  and  fSu  depend  on  a  and  6a 
only  in  a  bounded  subset.  Thus  we  shall  throughout  assume  that  a  and  8a  have  fixed, 
bounded  support. 

Let  D  C  IR*  be  open  and  bounded,  7  C  T*(D).  A  constant  C  is  said  to  depend  on  the 
Hs  D  HTmi( 7)-norm  of  w  G  ^'(IR1)  if  for  any  conic  neighborhood  T  of  7  there  exists  a  ip.d.o. 
Q  of  order  zero  with  ES(Q )  C  F  and  q  =  1  on  7  fl  {(a:,^)  :  |^|  >  1}  such  that  C  can  be 
bounded  in  terms  of  ||w||s,n  +  ||(5u>||r,n. 

The  following  is  the  statement  of  our  main  result  which  will  be  proved  in  the  sections 
which  follow. 

Theorem  1.1  Assume  either  that 

(i)  n  >  3,  5  >  max { 3  -f  n/2,  /  -f  n  —  1},  and  ao  G  Hs(JRn)  or 

(ii)  n  =  2,  s  >  max  {A,  l  +  1 },  and  a0  G  Hs  Fl  Hl+fl2({(x,  f)  :  8,n  =  0}). 

Then  under  the  assumption  (A),  the  following  estimate  holds 

\\DF(a0)8a\\l<C\\8a\\l+2fl  (1.8) 

where  the  constant  C  depends  on  the  Hs -norm  (n  >  3)  or  the  Hs  Pi  tf)  :  fn  —  0})- 

norm  of  u0  (n  =  2)  but  is  independent  of  8a. 


(1.5) 

(1.6) 


The  study  of  the  forward  map  is  motivated  by  the  inverse  pi'oblem  which  arises  in  reflec¬ 
tion  seismology,  oil  exploration,  ground-penetrating  radar,  etc.  Mathematically,  the  inverse 
problem  is  to  determine  the  coefficient  a  by  knowing  additional  boundary  value  conditions 
of  u.  Since  the  inverse  problem  is  just  to  invert  the  functional  relation  F,  we  are  naturally 
interested  in  all  the  properties  of  this  forward  map. 

To  understand  the  problem,  let  us  look  at  a  simple  exploration  seismology  experiment: 
Near  the  surface  of  the  earth,  a  seismic  source  is  fired  at  some  point  (point  source).  The 
seismic  waves  propagate  into  the  earth.  Since  the  earth’s  structure  varies  (as  do  its  physical 
properties)  part  of  the  energy  of  the  wave  will  be  reflected  back  to  the  surface  and  can  be 
measured.  The  inverse  problem  is  to  deduce  the  interior  properties  of  the  earth  from  the 
recorded  data. 

A  simple  model  of  this  reflection  seismic  inverse  problem  in  this  context  is:  given  data 
Fdata(%' ,t),  And  a  coefficient  <r(x)  so  that 


F(a)  =  Fdata 


or  perhaps  minimizing  the  error  (Fdata  —  F(a))  in  some  norm. 

Numerical  solution  of  this  problem  by  means  of  Newton’s  method  and  its  relatives  requires 
a  choice  of  Banach  space  structure  in  the  space  of  models  a  and  in  the  space  of  data  F(a) 
(see  e.g.  Kantorovich  and  Akilov  [15]),  in  such  a  way  that  F  is  singular.  The  simplest 
regularity  property  of  F  is  boundedness  of  DF ,  which  is  discussed  in  this  paper.  We  believe 
that  similar  arguments  will  establish  smoothness  of  F  and  allow  investigation  of  coercive 
properties  of  DF,  as  is  required  by  the  theory  of  optimization. 

The  most  efficient  relatives  of  Newton’s  method,  the  quasi-Newton,  conjugate  gradient, 
and  variable  metric  methods,  all  require  a  Iiilbert  space  structure  on  the  space  and  models. 
This  fact  accounts  for  our  reliance  on  the  Z2-based  Sobolev  spaces  in  this  work. 

When  the  spatial  dimension  is  one  or  c  and  a  depend  only  on  xn  (layered  problem)  there 
is  a  large  literature  available.  For  a  similar  problem  in  which  the  medium  was  assumed  to  be 
excited  by  an  impulsive  load  on  the  surface  {;r„  =  0}  instead  of  point  sources,  the  properties 
of  the  forward  map  have  been  studied  fairly  satisfactorily  by  Symes  and  others  (see  Symes 
[25]  for  references).  It  was  shown  by  Symes  that,  for  the  constant  wave  speed  case,  the 
forward  map  defines  a  Cl  —  diffeomorphism  between  open  sets  in  certain  Hilbert  spaces  by 
applying  the  method  of  geometrical  optics  together  with  energy  estimates. 

When  the  spatial  dimension  n  >  1  and  c,  a  depend  on  all  space  variables  (nonlayered 
problem),  very  little  is  known  in  mathematics.  Symes  [23,  24],  Sacks  and  Symes  [21],  Rakesh 
[IS],  and  Sun  [22]  have  some  partial  results.  The  difficulties  are  essentially  due  to  the 
ill-posed  nature  of  the  timelike  hyperbolic  Cauchy  problem  and  the  presence  of  nonsmooth 
coefficients.  For  the  one  dimensional  wave  equation,  both  coordinate  directions  are  spacelike, 
which  indicates  that  the  problem  is  hyperbolic  with  respect  to  both  directions.  Apparently, 
this  is  not  the  case  when  the  spatial  dimension  is  larger  than  one. 

Rakesh  in  [18]  looked  at  a  related  linearized  velocity  inversion  problem  with  constant 
density  and  point  sources.  Assuming  smooth  background  velocity,  he  obtained  some  results 
on  both  upper  and  lower  bounds  for  the  linearized  forward  map.  The  essential  observation  in 
Rakesh’s  work  is  that  DF  is  a  Fourier  integral  operator  (see  also  Beylkin  [7]).  Unfortunately, 
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the  calculus  of  Fourier  integral  operators  employed  in  Rakesh’s  work  is  not  applicable  to  the 
nonsmooth  reference  velocity  case  since  the  linearized  forward  map  is  a  Fourier  integral 
operator  only  when  the  reference  velocity  is  smooth. 

In  [23],  Symes  gave  a  pair  of  examples,  based  on  the  geometric  optics  construction,  which 
show  that  both  DF(1)  and  DF{ 1)_1  are  unbounded  for  a  slightly  different  problem.  As  the 
examples  show,  within  the  Sobolev  scales  no  strengthening  or  weakening  of  topologies  of  the 
domain  and  range  can  make  both  DF  and  DF~ 1  bounded.  This  fact  also  implies  a  strategy 
of  regularization:  Change  the  topology  in  the  domain  so  that  DF  becomes  bounded,  then 
ask  for  optimal  regularization  of  jDF-1  in  the  sense  of  best  possible  lower  bound  estimate 
for  DF.  In  both  examples  of  Symes,  the  unboundedness  was  caused  by  rapid  oscillation  of 
a  in  the  ^'-direction  or  the  tangential  directions,  hence  the  problem  is  actually  “partially 
well-posed”,  i.e.,  only  more  smoothness  of  the  coefficients  in  tangential  directions  (essentially 
grazing  ray  directions)  will  be  required  to  cure  the  difficulty.  This  might  be  the  main  reason 
the  anisotropic  Sobolev  spaces  IT"' "s ( ]Rn )  or  Hormander  spaces,  were  introduced  in  [21]  and 
[22]. 

In  Theorem  4.1  of  [21]  Sacks  and  Symes  showed  by  using  the  method  of  sideways  energy 
estimates  that  for  a  linearized  density  determination  problem  with  constant  velocity  and 
plane  wave  sources,  DF  is  bounded  from  H1,1  to  // 1 ,  provided  the  reference  coefficient  is  in 
H1,s  for  some  s  >  n  +  2.  They  also  proved  the  injectivity  of  DF.  Our  techniques  and  results 
are  quite  different  from  theirs.  We  intend  to  assure  the  optimal  regularity  of  the  timelike 
trace  under  different  hypotheses,  which  are  in  some  ways  weaker. 

There  remains  an  extremely  important  issue  to  be  addressed,  namely, 

What  is  an  appropriate  space  for  the  domain  of  DF  ? 

In  1983,  Symes  suggested  that  microlocal  restrictions  on  the  coefficients  might  regularize 
the  inverse  problem  (see  [24]).  In  some  sense,  this  was  confirmed  by  Bao  and  Symes  [2] 
where  we  were  able  to  prove  a  trace  theorem  for  the  solutions  of  general  linear  p.d.e.  with 
smooth  coefficients.  Roughly  speaking,  our  theorem  asserts  that  the  solution  will  belong 
to  Hs  along  a  codimension  one  hypersurface  if  it  belongs  to  Hs  in  a  neighborhood  of  the 
hypersurface  and  to  Hs+1  microlocallv  in  those  directions  where  the  p.d.e.  is  not  microlocally 
strictly  hyperbolic.  Note  that  we  gained  back  the  half  derivative  from  the  standard  trace 
theorem.  In  a  recent  paper  [3],  we  proved  a  similar  time  like  trace  regularity  result  for 
a  second  order  hyperbolic  equation  with  nonsmooth  coefficients.  It  is  obvious  that  the 
presence  of  nonsmooth  coefficients  will  introduce  new  singularities  to  the  solutions  so  that 
only  limited  initial  regularity  can  be  propagated.  A  crucial  step  in  [3]  was  to  develop  an 
extended  Beals-Reed  theorem  (Theorem  1  in  [6])  on  propagation  of  singularities. 

The  main  result  of  this  paper  is  a  boundedness  theorem  for  the  linearized  forward  map 
DF(cr0)  for  the  (sufficiently  regular)  nonsmooth  a0-  The  main  ingredients  of  our  proof  are 
the  method  of  energy  estimates,  a  microlocal  regularity  study  of  the  fundamental  solution, 
results  on  propagation  of  singularities,  several  trace  regularity  results,  and  a  useful  duality 
technique. 

It  is  known  that  in  their  applications  to  nonlinear  wave  equations,  most  of  the  results 
based  on  Rauch’s  lemma  (or  the  method  of  Fourier  analysis)  are  limited  to  relatively  weak 
singularities.  This  work  exhibits  that  to  some  extent,  strong  singularities  appearing  in  the 
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linear  wave  equation  ( e.g .  the  fundamental  solution)  can  also  be  tackled  by  this  Fourier 
analysis  method  with  the  help  of  a  duality  argument  and  the  progressing  wave  expansion. 
The  relation  between  the  coefficients  and  solution  with  strong  singularities  remains  to  be 
fully  understood,  especially  when  the  coefficients  are  less  regular. 

The  plan  of  this  paper  is  as  follows.  In  Section  2,  a  regularity  theorem  for  the  solution 
of  the  model  problem  is  established.  To  serve  this  purpose,  a  dual  problem  is  introduced.  A 
crucial  step  is  the  derivation  of  an  explicit  estimate  from  Homander’s  result  on  propagation 
of  singularities.  Section  3  is  devoted  to  the  proof  of  our  main  theorem.  The  main  ingredients 
of  our  proof  are:  microloca.1  regularity  analysis  of  the  solutions  for  the  transport  equations; 
a  microlocal  version  of  the  classical  trace  theorem;  and  regularity  study  of  the  dual  problem. 

Notation.  Throughout  this  paper,  the  reader  is  assumed  to  be  familiar  with  the  basic 
calculus  of  Pseudodifferential  Operators  (“  f.d.o.  ”)  as  stated  in  Taylor  [26]  or  Nirenberg 
[17].  A  classical  i 'p.d.o.  P  of  order  m  is  denoted  as  P  6  OPSm  with  its  symbol  p  €  Sm. 
ES(P)  stands  for  the  essential  support  of  operator  P.  W F(u)  denotes  the  wave  front  set 
of  a  distribution  u.  Hs  is  the  standard  L2-type  Sobolev  space  and  Hfoc  means  a  local 
Sobolev  space.  {£)  means  (1  +  |£|2)R2.  For  a  nice  discussion  on  microlocal  Sobolev  spaces 
Hs  Pi  IPmt(x o,£o),  we  refer  the  reader  to  Beals  [5],  see  also  Rauch  [19].  For  simplicity,  C 
serves  as  a  generalized  positive  constant  the  precise  value  of  which  is  not  needed. 

Warning.  When  the  reference  density  a0  is  smooth,  most  of  the  regularity  results  for 
the  forward  map  in  this  work  will  follow  more  easily  from  the  calculus  of  Fourier  Integral 
Operators.  For  a  standard  text  on  F.  I.  0.  we  refer  to  Duistermaat  [11]  or  Hormander  [13]. 
However,  this  technique  fails  with  the  appearance  of  the  nonsmooth  reference  density,  an 
assumption  important  in  this  work. 

2  Regularity  of  the  Fundamental  Solution 

Since  the  excess  pressure  u  in  the  model  equation  is  in  fact  the  fundamental  solution,  in  order 
to  study  the  regularity  of  the  forward  map,  the  regularity  of  the  fundamental  solution  must 
be  understood.  It  is  evident  that  the  real  obstacle  here  is  the  singular  right-hand  side  so 
that  none  of  the  propagation  of  singularity  results  could  be  applied  directly.  A  natural  way 
to  cure  this  difficulty  is  by  employing  the  Hadamard  theory  of  progressing  wave  expansion. 
We  refer  the  reader  to  Courant  and  Hilbert  [10]  or  Friedlander  [12]  for  a  detail  study  on 
the  method  of  progressing  wave  expansions.  According  to  Hadamard’s  construction,  the 
fundamental  solution  may  be  represented  as  a  sum  of  the  principal  part  and  remainder.  One 
can  then  study  the  remainder  by  the  Beals-ReecI  type  propagation  of  singularity  theorem. 
However,  a  great  drawback  of  this  idea  is  that  additional  regularity  is  needed  to  regularize 
the  remainder  term.  In  this  section,  we  develop  a  new  approach  based  on  the  method 
of  microlocal  analysis.  A  dual  problem  is  introduced  so  that  the  regularity  study  of  the 
fundamental  solution  may  be  transformed  into  regularity  study  of  the  dual  problem  which 
has  a  smooth  right  hand  side.  In  this  process,  a  crucial  step  is  to  derive  an  estimate  out 
of  Hormander’s  theorem  on  propagation  of  singularities.  As  one  might  expect,  with  the 
presence  of  nonsmooth  coefficients,  the  Rauch-type  results  and  some  commutator  results 
will  be  demanded. 
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2.1  Statement  of  result,  preliminaries 

Consider  a.  problem  obtained  by  integrating  the  model  problem  in  the  time  variable, 

(□  -  V<r0  •  V)n0  =  ,  (x,t)  G  IRn+1  f0  . , 

v0  =  0  ,  t  <  0  .  [  ’ 

The  following  is  a  regularity  theorem  for  the  fundamental  solution. 

Theorem  2.1  Suppose  that  1  +  n/2  <  s  and  gq  G  Hs(JRn).  Then  for  Vo  solving  the  equation 
(2.1),  l  <  s  —  n/2 

»o  €  ffL(U) 

where  U  —  {IRre  x  ( 0 ,  T i ) }  D  {/,  >  |x|}  (T\  >  0).  And  for  <j>  G  C^°(U),  the  following  estimate 
holds 

\\<fro\\i<C  (2.2) 

where  the  constant  C  depends  on  <b  and  1 1 rr0 j  | s  - 

In  order  to  establish  Theorem  2.1,  we  need  the  following  results.  The  first  was  originally 
established  by  Bony  [8]  and  was  extended  bv  Meyer  [16].  See  also  Beals  [4]  for  a  different 
proof. 

Proposition  2.1  Suppose  that  for  some  (.t0,£o)  G  T*(IR”)\0,  u  G  H*  Pi  Hrmt(x o,£o),  u/2  < 
s  <  r  <  2s  —  n/2,  and  g  G  C°° ,  then 

G  Hs  H  HTmi{x(h  to)  ■ 


We  also  need  a  Gardi rig's  type  inequality  concerning  the  microlocal  ellipticity. 

Lemma  2.1  Assume  that  Q\  G  OPSm 1 ,  Q2  €  OPS™2,  with  mi, m2  €  IR.  Furthermore 
assume  that  Q2  is  elliptic  on  ES{Q\).  Then  jor  any  r  G  IR,  S2  and  O'  two  open  bounded  sets 
o/IR"  with  fi  CC  W,  and  u  G  Q]°(0), 

IIQi«lkn  <C||Q2«||  s+mj  —m2  +  G||u||r,m. 

Proof.  Let  fR  and  fl2  be  open  sets  with  f!  CC  fli  CC  fl2  CC  ft'.  Construct  a  cut-off 
function  cf>  G  C£°(Q.'),  (p  =  1  on  fR,  and  <j>  =  0  on  f]'\ft2. 

The  assumption  Q2  is  elliptic  on  ES(Q  1)  implies  that  a  ip.d.o  R ,  a  parametrix  of  Q2  on 
ES(Q  1),  may  be  found  such  that 


Qi  IIQ2  —  Q 1  +  A 


(2.3) 


with  K  a  smoothing  operator. 

Having  defined  <j>,  we  can  now  rewrite 


QiRQ2u  =  Q\R(pQ2u  +  Q\R(l  -  <f>)Q2u  -v 
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It  follows  that,  for  any  r, 

\\Q1RQ2u\U  <  \\QiR<PQ2u\\s,n  +  \\QiR(i  -  <t>)Q2u\\sv 

<  C\\<pQ2u\\  s+mi  —  m2  ,!Rn  C||(l  -  (j))Q2u\\r<Q' 

<  C\\Q2u\\  s+mj- m2, fi'  "f”  ^  IMU*  • 

On  the  other  hand,  from  (2.3),  it  is  obvious  that 

\\QiRQ2u\\s,q  >  ||<2iu||s,n  -  ||Arn||s,n  . 

A  combination  of  the  above  discussions  will  complete  the  proof.  □ 

Finally,  the  proof  of  Theorem  2.1  demands  the  use  of  the  following  two  lemmas.  Lemma 
2.2  gives  an  estimate  based  on  Nirenberg’s  proof  [17]  of  Hormander’s  theorem  which  describes 
the  propagation  of  regularity  along  bicharacteristics.  With  nonsmooth  coefficients,  only  a 
limited  amount  of  regularity  propagates.  It  indicates  that  an  estimate  may  be  derived  near 
any  bicharacteristic,  hence  near  the  characteristic  variety  of  operator  □  =  d/  —  A.  We  then 
proceed  in  Lemma  2.3  to  argue  that  in  the  elliptic  region  of  the  operator  □  an  estimate  may 
also  be  formed.  With  a  concern  about  the  nonsmooth  <7o,  it  should  not  be  surprising  that 
both  results  require  a  commutator  argument. 

Let  IT  :  T*(Q,q)  — >  fio  denote  the  projection  of  T*(fi0)  onto  its  base  space. 

Lemma  2.2  Assume  that  s  >  n/2,  a0  €  //s+1(IRn).  Suppose  that  f3  is  a  null  bicharacteristic 
°fa>  (£(n£o)  €  0,  w  is  smooth  in  a  neighborhood  of  x0,  and 

Ow  -  Vcio-Vw  €  L2(l Rn+1) 

vanishes  in  a  neighborhood  of  11/3.  Then  there  exists  a  if.d.o.  B  of  order  zero  (essentially 
supported  near  j3)  so  that  for  any  <px  £  C,f’(IRn+1)  and  k  <  s  -  n/2  +  2,  there  is  C  >  0, 

||^i^lk<C||nu;—  Vdo-Vtullc 

Here  C  depends  on  crQ,  k,  B,  and  <j>lt  but  not  on  w. 

Proof  Let  fi  be  an  open  bounded  set  containing  supp(<f>  1)  and 

□  u»  —  V<T0  •  Vw  =  f  . 

According  to  Nirenberg’s  construction,  one  can  find  a  if.d.o.  Bq  of  order  zero  with 

(1)  b0  supported  in  a  small  conic  neighborhood  of  /?,  B0  elliptic  near  /?, 

(2)  El supp(bo)  fl  supp(f)  =  0,  and 

(3)  [O,  Bo]  6  OPS0. 
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Since  w  solves  (2.1),  the  method  of  energy  estimates  yields 


IMIi  <  c 


where  C  is  a  constant  depending  on  j j rr0 1 1 5  for  J  >  nf 2. 
Observe  that  from  (2.1), 


O[B0w  =  [□,  B0]w  +  [ B0 ,  Vtr0  •  V]w>  +  B0f  . 


Since  Ylsupp(bo)  D  stipp(  f)  —  0,  we  have 


B0f  =  0  . 


Now  energy  estimates  give 

H-^oHkn  £  (-'{ ||[a^  J30]iw||i,n  +  ||[-So;  Vco  •  V]u;|kQ)  .  (2-4) 

Since  [□,  Bo]  is  of  order  0, 

H[°,fioHI,.n<C’IMI,  <  C'll/llo  • 

The  third  term  in  (2.4)  may  be  estimated  by  applying  the  generalized  commutator  lemma, 
Lemma.  2.4  in  [3]  and  the  corresponding  estimate.  In  fact,  let  1  +  n/ 2  <  s0?  we  then  have 

||[5o,V(To-v«>]||lin<qH|1<c||/||o, 

where  C  depends  on  1 1 V cr0 1 1 s0  • 

Thus 

H-^oHkn  <  Colli  Ho  ,  (2-5) 

with  C0  depending  on  ||W0||4.0. 

Applying  Nirenberg’s  construction  once  again,  we  can  find  a  tp.d.o.  B 1  of  order  zero 
such  that  ES(Bi)  C  ES(Bo)( strictly),  B\  also  has  properties  (1)  and  (2)  above;  moreover 
€  OP S~x  and  B0  is  elliptic  near  ES(Bi).  From  (2.1)  and  B\f  =  0, 

B\W  —  [□,  B\\w  +  [Bi,  Vc0  ■  V]w  . 

If  we  write  down  the  energy  estimates,  alter  a  simple  ?/>.d.o.  cut-off  on  B\ ,  we  will  find 

\\BM\ln  <  C|HI?  +  C\\A\[Bi,  V<r0  •  VHknll#Hkn  , 
where  A\  €  OPS'0,  ES(Bi)  C  ES(Ai)  C  ES(B0),  B0  is  elliptic  on  ES(Ai),  and  at  =  1  on 

ES(Bi)n{(x,oAt  l>i}- 

Now  since  w  €  H 1  D  ES(B0)),  Lemma  2.4  in  [3]  again  implies  that  [B\,  V<7o  •  V]u;  € 
H1  n  H^f(ES(Ai))  and 

||Aj[Si,  Vcr0  •  V]u;||2,n  <  C(|H|i  +  H^iHkfi)  • 

Here  C  depends  on  1 1 Vcr0| |Sl  for  2  +  nj 2  <  si. 
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Because  of  our  construction,  Bq  is  elliptic  on  ES(Ai)]  therefore  Garding’s  type  inequality 
Lemma  2.1  leads  to,  for  any  real  r  and  0.  CC  hi 

IMHkn  <  C\\B0w\\2^  +C|Hlr  <  C\\f\\o 

by  (2.5). 

Therefore  we  have  shown  that 


I \BMkn  <C,||/||0, 

where  C\  depends  on  ||V<x0||Sl. 

We  can  continue  this  process  by  constructing  a  sequence  of  tf.d.o.  Bt  and  A{  (i  =  1,-  • 
•,  k  —  2),  such  that 

•  Bi  has  properties  (1),  (2),  and  [□,  B/\  £  OPS~\ 

•  ES(B^i)  C  ES(Ai-i)  C  ES(Bj ),  and 

•  Bi  is  elliptic  on  ES(Aj-i),  «;_]  =  1  on  ES{Bt^1)  Pi  {(x,£),  |£|  >  1}, 

•  Also 

||^i«!||i+2,n  <  Ct||/||o  , 

where  C;  depends  on  ||Vcr0|L  for  i  +  n/2  <  st. 

Eventually  we  conclude  by  choosing  B  =  Bk~2  so  that,  for  k  —  2  +  n/2  <  s, 

||-£Mkn  <  Gll/Ho 

with  C  depending  on  ||Va0||.,.  □ 

The  proof  of  Theorem  2.1  requires  a  slightly  different  form  of  Lemma  2.2. 

Corollary  2.1  Assume  that  s  >  n/2.  a0  6  Hs+1  (IR") .  Suppose  that  7  is  a  set  of  null 
bicharacteristics  0}  □,  (xo,£o)  €  7,  w  is  smooth  in  a  neighborhood  of  x 0,  and 

□  u-  -  V<t0  •  Yw  6  L2{ R,l+1) 

vanishes  in  a  neighborhood  oj  LI7.  Then  there  exists  a  if.d.o.  Q  of  order  zero  (essentially 
supported  near  7)  so  that  for  any  <7  £  C£°(IRn+1)  and  k  <  s  -  n/2  +  2,  there  is  C  >  0, 

\\<j>iQw\\k  <  C||nu’  -  Va0  •  Yw\\0  . 

Here  C  depends  on  a0,  k,  Q ,  and  <57.  but  not  on  w. 

Proof  For  every  null  bicharacteristic  of  the  set  7,  Lemma  2.2  indicates  that  a  if.d.o.  B  of 
order  zero  may  be  found  so  that 

WfaBwWk  <  C||Du;-  Vff0-  Yw\\0  . 

Now  Q  may  be  constructed  as  Q  =  B.  Moreover,  the  local  compactness  of  the  unit  sphere 
ensures  that  the  summation  is  finite.  □ 
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Lemma  2.3  Assume  that  s  >  n/2,  a0  €  Hs+1  ( 1R") .  Suppose  that  P  is  a  xf.d.o.  of  order 
zero  such  that  a  conic  neighborhood  of  its  essential  support  is  contained  in  the  elliptic  region 
of  the  wave  operator  □.  Assume  also  that 


□?(<  -  V«70  •  Vw  €  X2(IRn+1) 

vanishes  in  a  neighborhood  of  Tip.  Then  for  any  fi  6  Qf  (IRn+1)  and  k  <  s  —  n/2  +  2,  there 
is  C  >  0  so  that 

H^iPn’Ik  <  C\\Ow  -  V(J0  •  Vio||0  • 

Here  the  constant  C  depends  on  cq,  k.  P,  and  f>\,  but  not  on  w. 

Proof.  The  proof  is  based  on  the  same  type  of  bootstrap  arguments  as  in  the  proof  of  last 
lemma. 

Let 

□  »’  —  V<To  •  Vto  =  /  .  (2-6) 

Assume  that  ft  D  supp(<f> i)  be  a  bounded  open  set.  From  the  support  assumption  on  p,  we 
see  that  P f  =  0.  Hence,  by  applying  P  to  both  sides  of  (2.6),  we  find 

OPw  =  [□,  P]w  +  [P,  V<t0  •  V]tc  +  V<j0  ■  VPw  .  (2-7) 

Now  since  □  is  elliptic  in  a  small  conic  neighborhood  of  ES(P),  there  exists  a  ijj.d.o.  Pq  of 
order  zero,  such  that  ES(P)  C  ES( Pq),  Pq  is  elliptic  near  ES(P ),  and  □  is  elliptic  in  a  small 
conic  neighborhood  of  ES(Pq).  From  the  ellipticity  of  P0D  on  ES(P),  Lemma  2.1  gives,  for 
any  real  number  r  and  fi  CC  fi', 

ll-Ptelk.n  <  C\\P0nPw\\k-2,si'  +  C|Mk  , 

or  from  (2.7) 

Ili’Hkn  <  C{\\P0p,P]w\\k-w  +  ||Po[P,  V<70  •  V]to|U_2,n,  +  ||P0V<70  ■  VPu>|U-2,0')  . 
Therefore  an  application  of  Lemma  2.4  and  the  generalized  Rauch’s  lemma  in  [3]  yields 

ll^ielk.n  <  C  i  ||Po'U’||fc_i,n'  +  P b  ( 1 1  w  1 1 1  +  ||Pow|U-2,fl') 

+(M\H\i  +  ||Po^lU-i,n') 

—  P'lL/llo  +  C\\P0w\\k-l,n'  ■ 

Here  constants  C2  and  P3  depend  on  ||V<70||S  for  k  —  2  +  n/2  <  s. 

Thus  the  bootstrap  arguments  on  Pq  will  accomplish  the  proof.  □ 
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2.2  Proof  of  Theorem  2.1 


We  study  the  regularity  of  v0  through  its  dual  problem.  To  simplify  the  arguments  on  the 
dual  problem,  we  make  use  of  the  symmetric  form  of  (2.1)  by  introducing  p(x )  =  e_<r°.  Then 
(2.1)  becomes 

□,.-o  =  [by  -  V  •  (iv)].,„  =  -S-^{t)S(x) 

p  p  p  (2.8) 

Vo  =  0  t  <  0  . 

Now  let  us  look  at  a  dual  problem  to  (2.8), 


□  lU7  =  \-d?  —  V  •  ( -V)]n>  =  T 

P  P 

(2.9) 

w  —  0  t  »  T]  , 

where  T  €  Qf 

5(fi)  with  n  =  { ]R."  x  (o.T\)}  n  {t  >  |x| 

+  Co},  for  e0  >  0  small. 

Note  that, 

this  equation  may  be  reformulated  as 

□  ]  W  =  Uw  —  V (To  ■  Yw  = 

e_<T°T 

(2.10) 

tv  =  0  t  »  T\  . 

Thus  if  we  can 

show  that  for  any  T  G  6'0x,(O) 

|(^-oCk)|<C||vk||0 

) 

(2.11) 

then  it  can  be 

concluded  that 

II^WIIO.Q  <  C  . 

(2.12) 

From  (2.8), 

integration  by  parts  leads 

1(^-0.^  )|  =  u,)l 

=  |(  !«(<)«(*),  # 
n 

-(n— 1)/2  vi 
'  ic)| 

<  C|(^-(n-1)/2u;)(0,0)|  . 

The  trace  theorem  (see  for  example  [26])  yields  that 

|(5tV®)|<C||^HI/+1  (2.13) 

with  <f> i  €  Co°(fii),  Hi  a  small  neighborhood  of  the  origin  and  fii  fl  supp(ty)  =  0. 

Construct  two  ib.d.o.  Q\,  Qi  G  (9/CS'°(IR”+1),  such  that 

•  Qi  +  Q2  —  R',  where  R  is  an  elliptic  d’.d.o.  of  order  zero  in  fh; 

•  Hsupp(qi)  fl  supp(  T )  =  0,  for  i  —  1,2; 

•  ES(Q 2)  is  a  small  conic  neighborhood  of  set  of  null  bicharacteristics  of  the  wave  oper¬ 
ator  □  passing  over  f C; 
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•  Q i  is  microlocally  smoothing  on  the  null  bicharacteristics  passing  over  fli . 

Therefore,  with  (2.13),  we  have 

|(0tV*)|  —  +  C\\Q2(j)2w\\i-\-i,n1  (2-14) 

here  the  expression  makes  sense  because  the  domain  of  dependence  for  w  and  the  pseudo-local 
properties  of  Q\  and  Q 2. 

Now,  we  can  apply  Corollary  2.1  to  obtain  that 

<C||*Ho.  (2.15) 

Lemma.  2.3  yields 

IIQ2«'II/+i,o1  <<?ll*llo  (2.16) 

where  the  constants  here  depend  on  |  i cr0 1 1 s - 
Therefore,  we  have  shown 

|(c>tw^)|  <  C||^1«7||,+1  <  CII^H  .  (2.17) 


To  complete  the  proof,  we  have  one  more  step  to  go,  that  is,  to  estimate  all  the  x- 
derivatives  as  well  as  the  corresponding  mixed  derivatives  of  vq  up  to  order  of  l.  This  may 
be  done  by  introducing  another  dual  problem.  Here,  we  demonstrate  this  method  by  looking 
at  d‘x  Vo  (i  =  1,  •  •  *,n).  Tlie  rest  terms  can  be  estimated  by  the  same  fashion.  It  suffices  to 
show  that  for  any  T  €  C£°(f!) 

|(4(Co.^)|  <C||T||0  (2.18) 

for  i  =  1,  •  •  •,  n. 

Let  iv  solve  the  dual  problem  (2.9),  u>,  solve  the  following  problem 

□  iu;t  =  {-d'f  -  V  •  (-V)]tei  =  [dl  Diju;  (  , 

P  P  (2.19) 

iv  —  0  /  >>  I\  . 


We  have 


l(3S,*)|  =  |(u0,^nito)| 

<  |(niUo,^, .«•')!  +  |(OiV0,W7i)| 

<  C|(^r9r("",)/2rc)(0,0)|  +  C|(at-(ra-1)/V)(0,0)|  . 

Again,  the  trace  theorem  gives 

K^UoAlOl  <  C||^iu>||/+1  +  CH^itWilli 

<  C 1 1  <pi  te|  |/-)-i 

where  in  order  to  obtain  the  second  inequality,  we  have  applied  the  standard  energy  estimates 
and  a  commutator  argument.  From  (2.17),  the  estimate  (2.18)  can  then  be  proved.  □ 
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3  Proof  of  the  Main  Theorem 


Our  goal  in  this  section  is  to  determine  the  appropriate  hypotheses  under  which  DF(aa), 
the  linearization  of  F  about  a  reference  state  <7o,  is  bounded.  Our  proof  is  based  on  the 
microlocal  regularity  analysis  of  the  transport  equations  and  regularity  study  of  the  dual 
problem. 

Recall  the  linearized  problem  corresponding  to  the  reference  state  (u0,cr0),  for  ( t,x )  € 
IRn+1,  x  =  (x',xn), 

( □  —  Vao  •  V)8u  —  V(5cr  •  Vu0  i , 

8u  —  0  ,  t  <  0  ,  [  } 

where  uq  is  the  solution  of  the  model  problem  corresponding  to  the  reference  density  ao.  The 
linearized  forward  map  can  be  defined  as 

DF(<jq)6ct  =  ( <j>8u )  |xn=o  ,  (3.2) 

where  <fi(x,t)  £  C(j°(IRn+1 )  is  supported  inside  the  conoid  {t  >  | } ,  and  near  {xn  =  0}. 

Once  again  we  consider  a  related  problem, 

(□  -  V(70  •  V)  v  -  VSa  ■  Vv0  ^  o\ 

v  =  0  ,  t  <  0  ,  ^  ‘  ’ 


where  8u  —  dt  2  v  and  v0  solves 


(□-V<T0  -V)ro  =  8-’-^-(t)6(x) 

no  =  0  ,  t  <  0  . 


Observe  that  for  /  €  IR., 


\\DF{a0)8a\\i  =  ||  (<f>8u)  U„=o||/ 

<  C||(H  l.n=o  Ih,  (3-5) 

where  l\  denotes  /  +  (n  —  l)/2.  Thus  the  real  challenge  here  is  to  get  an  appropriate  trace 
regularity  estimate  for  v  on  a  time-like  hypersurface  {xn  =  0}. 

Before  getting  into  the  details  of  the  proof,  let  us  first  make  the  following  general  remarks 
on  this  theorem: 

The  estimate  (1.8)  has  a  similar  form  to  a  Rakesh’s  theorem  (Theorem  2.5  in  [18]). 
Actually,  we  conjecture  that  a.  formal  extension  of  our  proof  here  could  lead  to  an  elementary 
proof  of  his  theorem.  On  the  contrary,  the  principal  tool  of  Rakesh’s  proof,  calculus  of  Fourier 
integral  operators,  is  not  available  when  the  reference  density  is  nonsmooth. 

Our  approach  enjoys  the  beauty  of  the  method  of  energy  estimates,  that  is,  it  possesses 
useful  information  on  various  parameters  involved  in  the  estimates. 
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3.1  Some  useful  results 

Hadamard’s  construction  leads  to  the  progressing  wave  expansion  for  n0, 

Co  =  bkSk :{t  -  t{x))  +  RVo(x,t)  (3.6) 

k=0 

where  t(x)  =  |.r|,  So  is  the  Heaviside  function,  S'k  =  Sk- 1  ( k  >  1),  and  RVo  vanishes  at 
t  =  t(x).  Moreover  { b *}  solve  the  transport  equations,  for  k  =  1,-  •  -,s, 

2Vr  •  Vb()  +  (At  +  Vr  ■  Va0)b0  =  0  (3.7) 

2Vr  •  Vbk  +  (At  +  Vr  •  Va0)bk  =  A6,_,  +  Va0  •  Vbk^  .  (3.8) 

From  (3.3),  Hadamard’s  construction  also  gives  the  progressing  wave  expansion  of  v, 

S 

v  =  ]T  akSk(t  -  t(x))  +  Rv(x,  t)  ,  (3.9) 

k=0 

where  r( x)  —  |.r|,  ,S'o  is  the  Heaviside  function,  S'k  =  Sk~u  Rv  vanishes  a,t  t  =  t(x),  and  {ak} 
solve  the  transport  equations,  for  k  —  0,  •  •  -,s  —  1, 

2 Vr  •  V«o  +  (At  +  Vr  ■  Va-0)a0  =  —boVr  ■  VSa  (3.10) 

2 Vr  •  Vak+ 1  +  (Ar  +  Vr  •  Vcr0W+i  =  A ak  +  Va0  ■  Vak  +  V^cr(V6fc  -  bk+1Vr)  .  (3.11) 

In  fact,  the  equations  (3.10)  and  (3.11)  are  the  first  order  perturbation  of  the  equations  (3.7) 
and  (3.8),  respectively. 

For  convenience,  we  introduce  a  function  q  =  ao/2  +  qo  with  Vr- V</o  =  At/2.  Thus  away 
from  the  origin,  q  is  nothing  more  than  a.  smooth  perturbation  of  ao/2.  Then  the  transport 
equations  (3.7),  (3.8)  may  be  transformed  to  equations 

Vr-V&oe"  =  0  (3.12) 

V r  •  Vbke'J  =  ( A6fc_a/2  +  W0  •  V6fc_x/2)e«  ,  (3.13) 

for  k  =  1,  •  •  •,  /i. 

Similarly,  equations  (3.10)  and  (3.11)  may  be  transformed  to  equations 

Vr  •  Va0e9  =  Vr  •  V(6aboeq)  (3.14) 

Vr  •  Va,.e"  =  (A«,..1/2  +  V<70  ■  Vafc.1/2)e?  +  Vfo(V5fc  -  &H1Vr)e?/2  ,  (3.15) 

for  k  =  1,  •  •  •,  li,  where  to  obtain  the  first  equation,  we  have  used  the  equation  Vr  •  Vfroe9  =  0. 
Observation.  The  right  hand  side  of  the  equation  (3.13)  may  be  rewritten  in  terms  of  bk~ieq 
as 

A(6,_ief0/2  -  I VqW^eyo  -  6fc_ie»A q/2  -  V(6fc_ie")  •  Vq0  +  h^e'Vq  ■  Vq0  .  (3.16) 

Similar  observation  may  be  made  for  the  right  hand  side  of  (3.15). 
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Observe  that  all  the  transport  equations  have  the  same  principal  part  Vr  •  V  which  is 
a  smooth  vector  field.  Therefore,  in  order  to  understand  the  regularity  of  the  solutions  to 
(3.12)-(3.15)  it  is  essential  to  study  the  properties  of  this  smooth  vector  field.  The  following 
is  a  microlocal  regularity  result  for  the  solutions  of  transport  equations.  It  indicates  that  a 
refined  regularity  result  can  be  achieved  microlocally  away  from  the  characteristic  variety  of 
the  wave  operator,  in  contrast  to  a.  direct  application  of  the  energy  estimates.  This  result 
will  be  used  frequently  throughout  the  proof  of  Theorem  1.1. 

Lemma  3.1  Let  V  be  a  smooth  vector  field.  Assume  that  u  is  a  smooth  function,  Q  £  OPS0 
supported  away  from  Char(V),  and  where  Si  is  a  bounded  open  set.  Then  there 

exist  Q'  £  OPS0  also  supported  many  from  Char(V )  and  ft,  f  £  C'ft (IT)  with  SI  CC  Si1, 
such  that  for  any  r  £  1R, 

||Q<MLn  <  C\\Q'<f)'Vu\\s-isi'  +  C'llHU  • 

Proof.  Let  Q\  be  another  fi.d.o.  of  order  zero,  such  that  Qi  is  elliptic  on  ES(Q)  and 
ES(Qi)  is  away  from  Char(V).  It  follows  from  Lemma  2.1  that  there  exists  Qi  DD  Si  an 
open  bounded  set  such  that 

||Q<HLn  <  CIlQiV'^ulls.n,  +  C||<Hkni  (3-17) 

<  T'||Qi[i’.  d]'(/||s-iQi  +  CHQi^Vulls-i.ni  +  •  (3.18) 

Now  choose  <f\  £  C£°(fii)  and  <j>x  =  1  on  supp(<f>),  then 

Q i [V, c6] =  Qi[V,<p\(piu  -  Q\fi\u 

where  Q\  =  Q\[V,  <j>]  is  a  if.d.o.  of  order  zero  and  ES(Q i)  is  away  from  Char{V). 

Thus 

IK,?<H|5,q  <  CHQkMIs-ia  +  CWQxfiVuWs-i^  +C||^t<||r,fi1  . 

We  can  then  repeat  the  above  process  to  get  a  similar  estimate  on  •  In  fact, 

one  can  construct:  Q2  £  OPS0,  Q 2  is  elliptic  on  ES(Qi)  and  ES(Q2 )  is  away  from  CharfiV)-, 
<t>2  £  Cff{Sl2),  (p 2  =  1  on  supp(fi),  ill  CC  Sl2;  Qi  =  Qi^fifi). 

Then 


IIQkMUa  T  C\\Q2<p2u\\s-2,n2  +  OWQifiV u\\s^2,u2  +  C|kin||r,n2  . 

Following  the  same  idea,  one  may  keep  going  to  get  similar  estimates:  (the  i-th  step) 
Construct:  Qi  £  OPS0,  Q ,  is  elliptic  on  ES(Qi-i)  and  ES(Qi)  is  away  from  CharfV ); 
fi  £  Cff(Sli),  ft  =  1  on  supp(6,-i )-  SI*-— i  CC  Sip,  Qi  =  Q,[V,  ft],  0  <  i  <  k  —  1.  (k  is  chosen 
to  satisfy  s  —  k  <  r). 

Then 


1 1  QifiiU  1 1  s_  t,fi;  T  O  1 1  (Ji+l  4>i+  1  n  |  |.s 


+  1 


+  C\\Qt+i4>iVu\\s 
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Combining  all  of  the  estimates  above,  we  have 


k-l  k- 1 

+  C  Y2  ||Q4+i<M/ulls-i-i,ft,+i2  T  C  ll^iwlkn,+i  • 

*=0  i— 0 

In  order  to  get  the  desired  estimate,  some  further  simplifications  are  needed.  Clearly,  for 

j  =  0,  •  •  •,  k  -  2, 

ll<Mkn;+i  <  C\\<f>j+iu\\rflj+2  , 

hence 

k-i 

<  C\\(f>k-iu\\r<Qk  . 

i=0 

Since  s  -  k  <  r,  ||(5fe^A.-«||5-it,nfc  <  C  \  \ <pkU  1 1 r,nk ,  therefore, 

k- 1 

IIQ<HU  <  qi<MU  +  J2  \\Qi+i<i)iVu\\s-i-i,n,+1  • 

i= 0 

Construct  a.  xji.d.o.  of  order  zero,  Q\  with  the  following  properties: 

Q'  is  elliptic  on  u£r01jE?S'(Qt-+i<p1-);  ES(Q')  is  away  from  Char(V).  Lemma  2.1  then  yields 
that 

\\Qi+i<f)i\  u | |s- 1 ,fj;.  =  | \Q  i+ififkV  |s-i,n^-  (3.19) 

<  CWQ'faVuW'-w  +  C\\<t>kVu\\r.hw  (3.20) 

where  f 7*  CC  Cl\  i  =  0, 1,  •  •  •,  k  —  1.  Therefore,  if  we  construct  <j>  G  C£°( O')  in  such  a  way 
that  ar<j)  >  0  on  supp(<f>k),  then 

ll^llr  — 1,0'  hi  C||<Mlkfi'  . 

Let  <j>'  be  <j>k,  we  finally  obtain 

[|Q<HLq  <  CWQ'&VuWs-i#'  +  C\\4>u\\r<Q> 

which  completes  our  proof.  □ 

We  also  need  the  following  standard  result  for  hyperbolic  p.c/.e,  as  well  as  the  estimates 
involving  in  its  proof.  See,  for  example,  Chazarain-Piriou  [9]  for  the  idea  of  the  proof.  The 
following  is  the  version  stated  in  Beals  [5]. 

Lemma  3.2  (Linear  Energy  Inequality)  Let  p(x,D )  be  a  partial  differential  operator  of  or¬ 
der  m  on  IRn+1 ,  strictly  hyperbolic  with  respect  to  the  plane  {xn+1  =  0},  and  let  u  satisfy 
p(x,D)u  =  f(x).  If  f  G  L//f“n,+1  (1R'!+1 )  and  u  G  Hfoc(x  :  |xn+i|  <  e)  for  some  e  >  0,  then 

ueHfj  rcG1). 

Next,  we  present  a  trace  regularity  result. 
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Lemma  3.3  Assume  that  s  >  3  +  »/ 2,  1  <  /j  <  s,  and  v  solves  problem  (3.3)  then  there  is 
a  4> o  €  C™  supported  near  supp(cp)  such  that  the  following  estimate  holds, 

||(oi>)U„=olk  <  C\\<j>0v\\h  ,  (3.21) 

where  C  is  a  constant,  depending  on  the  IIs  Pi  Hlffl  (K)-norm  of  gq,  but  is  independent  of 
8a . 

Proof.  This  lemma,  is  a  direct  application  of  Theorem  3.1  in  [3]  by  taking  into  account  of 
the  fact  that  cf>  and  8a  have  disjoint  supports.  P 


3.2  Microlocal  version  of  trace  theorem 

The  classical  trace  theorem  in  Sobolev  spaces  characterizes  the  regularity  of  a  distribution 
restricted  to  a  hypersurface.  Dealing  with  inverse  problems,  one  always  has  to  face  a  difficult 
but  crucial  question:  When  does  the  restriction  operator  commute  with  another  operator 
of  interest?  The  result  in  this  subsection  indicates  that  a  simple  microlocal  trace  theorem, 
which  not  only  works  on  the  space  restriction  but  also  on  the  phase  restriction  (be.  a  trace 
theorem  on  cotangent  bundles),  may  lead  to  a  way  to  cure  this  difficulty. 

Let  K  be  a  conic  set  in  lRn.  i  :  x  6  IR"  — >  (x,  0)  €  IRn+1.  Define  a  semi-norm:  for  7  a 
conic  set  of  IRfe  and  u  €  ( 1R k ) , 

M,.,  =  ( /  <(f|fi({)l2(02')1/2  • 

■tie- 

Then,  a  proof  of  the  classical  trace  theorem  (see  e.g.  in  Taylor  [26],  pages  20-21)  implies  the 
following  inequality. 

Proposition  3.1  For  s  >  1/2.  u  €  6'o°(lRn+1), 

\i*'M  |  K.s-l/2  <  C’|u|xx]R,s  • 

Thus  the  map  i~  may  be  extended  to  be  a  bounded  map  from  tlsmt{x  x  IR,  K  x  IR)  to 
Hme1/2(xiK)’  provided  s  >  1/2. 

Let  n2  be  the  projection  map  to  the  frequency  space  (or  the  second  factor).  We  may 
reformulate  this  result  in  terms  of  tp.d.o.. 

Proposition  3.2  If  P\  is  a  xp.d.o.  of  order  zero  in  IR",  with  II 2 -£',$'( Pi )  C  K,  then  there 
exists  a  if.d.o.  P2  of  order  zero  in  IR’I+J ,  and  Tl2  FS(  P2)  C  K  x  IR  ,  such  that  for  s  >  1/2, 
u  £  C~(fl)  with  Q  an  open  bounded  subset  of  IRU+1 , 

||Pi?'*'i<||s-i/2,n0  —  C||-f2w||s,n  , 

where  i*  again  denotes  a  restriction  operator  to  a  codimension  one  hypersurface  and  Qq  = 

i*n. 
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The  above  results  together  with  our  Garding’s  type  result  Lemma  2.1  yield  a  microlocal 
version  of  trace  theorem. 

Lemma  3.4  Assume  that  E  is  an  elliptic  operator  of  order  m  in  IRn+1  x  K  x  IR,  P  G 
OPS°(lRn)  and  U.2ES{P)  C  I\  .  Then  for  s  >  1/2,  u  G  Cq°(Q)  where  f l  and  ft'  are  open 
bounded  subsets  o/]Rn+l  with  Tl  CC  O',  and  =  i*Tl, 

||^:*«||,-,/2.Qo  <  C\\Eu\\t-m,0'  +  C||«||r,n. 


for  any  r  G  IR . 

A  result  on  propagation  of  singularities,  see  Proposition  1.3.3  in  Duistermaat  [11]  or 
Theorem  8.2.13  in  Hormander  [14],  demonstrates  the  relation  between  the  wavefront  of  the 
restriction  of  a  distribution  and  the  wavefront  set  of  its  own.  However  the  result  does  not 
directly  lead  to  any  explicit  bound.  Here,  in  this  direction,  we  present  a  result  together  with 
an  estimate.  The  proof  follows  directly  from  Lemma  2.2  and  Lemma  2.3. 

Lemma  3.5  Assume  that  s  >  n/2,  cr0  G  //s+1( lRn).  Suppose  that  w  is  smooth  near  {t  =  0} 
and  Dto  —  Vaa  ■  Vw  G  L2(IRn+1).  Then  there  exists  an  elliptic  xf.d.o.  B  of  order  zero ,  such 
that  ES(B)  is  contained  in  Cy,  a  "cylindrical”  conic  neighborhood  of 

{(•M,£,cr)  G  7”*IR"+1\0,  f2  -  | m [ 2  =  0,w  =  Vr  •  (} 

along  uj  direction ,  and  the  symbol  of  B.  b  satisfies 

Tlsupp(b)  fl  supp{ □«)  —  V(7o  •  Vrt>)  =  0  . 

Then,  for  any  <pi  G  C^°(Kn+1)  and  k  <  s  —  n/2  +  2, 

1 1 c6 1  /?».) 1 1 <  C'l 1 1  □  uj  -  V (7q  ■  Vu>||0  , 

where  the  constant  C  depends  on  <tq,  k,  B,  and  fix,  but  not  on  w. 

3.3  Proof  of  Theorem  1.1 

Recall  the  assumption  made  in  Sect  ion  1. 

(A)  supp(Sa)  C  {v„  >  <:}  , 
for  e  >  0  small. 

For  simplicity,  we  shall  also  assume  that  l\  is  an  integer.  Without  any  further  difficulty, 
the  proof  may  be  extended  formally  to  cover  the  general  case. 

Let  <p  G  C'fi  be  supported  inside  the  characteristic  surface  and  the  set  {xn  <  e/2}. 
Multiplying  <j>  to  both  sides  of  equation  (3.3),  we  have 

•  lor  =  dV<70  •  Vr  +  [□,  <f>]v 
r  =  0  ,  t  <  0  . 


(3.23) 


(3.22) 
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Here  we  ha.ve  used  the  fact  that  according  to  the  assumption  (A),  (j>  and  da  ha.ve  disjoint 
supports,  so  that  4>Y8a  ■  Vvo  =  0. 

Having  Lemma  3.3,  the  estimate  of  ||(<?!w)L,=olk  maY  be  reduced  to  estimating  ||^0u|k- 
For  the  sake  of  simplicity,  we  shall  only  discuss  the  estimates  of  the  t-derivatives  of  v.  The 
rest  may  be  estimated  as  in  the  proof  of  Theorem  2.1. 

The  rest  of  this  section  is  devoted  to  the  estimate  of  ||</>ou|k- 

We  study  the  regularity  of  v  through  its  dual  problem.  Once  again,  we  look  at  the 
symmetric  form,  for  p(x)  —  e~a° . 


□  m  -  \-df  -  V  •  (-V)K  -  -V6cr  •  Vuo 
P  P  P 

v  =  0  t  <  0  . 


(3.24) 


A  dual  problem  to  (3.24), 

a\w  =  [~r)f  -  V  •  (~Y)]w  =  $ 
w  =  0  t  »  T\  , 

where  T  €  with  Q  =  {1R'1  x  (0,T))}  fl  {f  >  |r|  +  e},  for  a  small  e  >  0. 

Thus  if  we  can  show  that  for  any  T  €  Qj°(h) 

\(d^v2,n<C\\Sa\\h  ||T||o, 

then  it  can  be  concluded  that 

||#t*||o,o<C||Hk- 

Green’s  identity  and  integration  by  parts  lead  to 

T)  =  (^‘c.Dju?) 

=  —  2  /  —Idi1  v  — — w  —  w—di'vjds 

Ji-r(x)  p  On  on 

Let  us  look  at  the  first  term.  From  equation  (3.3), 

( □i(d(1  v,  w)  —  (VScr  ■  Yd{‘v0,  w)  . 

Thus 

|(D1dlS\'U.’)l  <  C'll^ll  ||dNo||o,n  IMkn  . 

Applying  Theorem  2.1  and  energy  estimates  to  the  right  hand  side,  we  have 

|(D1dt,1t;,u;)|  <  C||HI  ll^ll 

with  the  constant  C  depending  on  ||cr0||s,  for  s  >  max{\  +  n/2,l\  +  nf 2}. 

Using  Hardamard’s  construction,  one  may  get: 

($t  ’^  )  |  i = r  ( a- )  )  {&t  'e)|t=r(.r)  *Ri+l,  ( )  |  f =-r(a;)  9xal\  T  ®q+l- 


(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 
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Denote 


h 


1  ()  , 

-Otlv  —  wds 
t=r(x)  p  on 


-L 


1  d 


t—r(x)  p  9n 


(d{'  v)wds  . 


Observation.  The  integrands  are  all  compactly  supported  even  though  neither  v  nor  w  is 
compactly  supported. 

Since  the  above  integrals  have  similar  forms,  we  will  only  give  a  detail  estimate  of  I\  by 
making  the  following  remark  on  the  estimate  of  I2. 


h  =  -2  /’  -Vr  •  V(dlfv\t=T)wds  +  [ 

Jt  =  T(x )  p  Jt=^ 


(x)  p 


dlt1+1vwds 


Now,  by  integration  by  parts,  the  first  term  of  the  right  hand  side  can  be  studied  just  like 
Ii]  while  the  second  term  may  also  be  handled  in  a  similar  fashion  as  for  I\. 


3.4  The  estimate  of  I\ 

A  function  <j>{x)  G  C'^°(0),  Q  an  open  bounded  set  in  IR/1  and  D  3  supp{d\v\t=T-§^w\t-T} •, 
may  be  found  because  the  integrand  of  / 1  is  compactly  supported,  so  that 

./,  =  /  (j)-dltv—wds  . 

■h=T(r)  p  on 

Now  construct  two  if.d.o.  Qj,  Q>  <G  OPS°{ lRn) 

•  Ql  +  Qi  =  P, 

•  ES(Q 2)  is  a  small  conic  neighborhood  of  {Vr  •  £  =  0}; 

•  Q2  s  symbol  q2  =  1  near  {Vr  ■  £  =  0}  ft  {(a;,f),  |£|  >  1}. 

Recall  the  definitions  of  Q 1  and  Q 2.  We  can  again  rewrite  I\  as,  for  a  <fti(x)  G  with 

ft  C  D', 

h  =  /  QiO-dl1  vQxEi-^-wds  +  /  Qif-d^vQ^i—wds  +  [  Q^-d'^v^—wds 

J  t—T  p  On  J  t—T  p  On  Jt=T  p  On 

=  I l,a  +  h,b  +  h.c  • 

Now,  we  estimate  / 1,„,  I i.j,,  and  I i,c  separately. 

(a)  The  estimate  of  /l  a 

Clearly,  in  order  to  estimate  it  suffices  to  estimate  ||  and  \\Qi<t>i-§^w\t=T\\. 

These  terms  may  be  handled  by  the  following  propositions. 

Proposition  3.3  The  following  cslimate  holds 

WQ\4>bke.q\\o,n  <  Ck 

where  ft  CC  Cl  and  Of  depends  on  j | rr0 1 1 1 - 
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Proof.  The  proof  is  based  on  an  application  of  Lemma  3.2.  We  also  have  to  use  variants  of 
Rauch’s  Lemma,  Schauder’s  Lemma,  together  with  Proposition  2.1. 

In  fact,  Lemma  2.2  leads  to,  for  any  r, 

llQi^Ww  <  CWQ'^Vt  •  V6fce’||_ll0l  +  C\\fabke%,ai 

where  Q  CC  -Qj  ,  <j>[,  <j>\  G  ),  and  Q'  supported  away  from  Char(Vr  ■  V).  From  (3.13), 

it  is  easy  to  represent  the  right  hand  side  in  terms  of  bk-\eq .  Therefore,  Rauch’s  Lemma  and 
Proposition  2.1  yield 

WQppb^Wosi  <  C\\Q\fabk-ie%fll  +  C||MtlU 


where  s  >  1  +  n/2  and  the  constant  C  depends  on  ||<r0||s. 

One  can  then  continue  this  process  by  applying  Lemma  2.2  to  estimate  WQ^fcbk-i^Wifii- 
After  /  similar  steps,  we  can  write  the  following  estimate,  for  any  r, 

k 

IKM^IIo.u  <  C\\Q'ld>,b0eq\\ku  +  £C’||M||r,fit  (3.30) 

i= 1 


where  C  depends  on  | |<r0 1 [s  for  s  >  I  +  n/2,  fl  CC  Oi  CC  •  •  •  CC  Ofc,  and  <j>i  G 
We  next  discuss  why  the  smoothing  terms  ||<A6(||r,Qt  may  be  negligible. 

Note  that  the  transport  equations  are  hyperbolic  along  the  A  direction  (A  =  |ar |) .  Actually, 
introducing  polar  coordinates,  we  then  get  Vr  •  V  =  Jy  (A  =  |x|).  Thus,  A  may  be  treated 
as  the  “time”  variable  for  a  standard  hyperbolic  problem.  The  right  hand  sides  of  transport 
equations  may  be  expressed  under  the  polar  coordinates  as 


dbkeq 

d\ 


(3.31) 


Because  of  the  main  assumption.  a(x)  is  constant  near  |,r|  =  0,  that  is,  v0  is  the  fundamen¬ 
tal  solution  for  the  same  equation  with  constant  coefficients.  The  Hadamard  construction 
yields  that  bk  (k  =  0,  •  •  •,/])  are  smooth  function  near  |.r|  =  0  (or  A  =  0),  see  [20]  or  [1]  for 
more  discussions.  Lemma  3.2  then  becomes  applicable. 

Eventually,  there  exists  a  constant  Ck  depending  on  ||cr0||fc,  so  that 


IKM^H  <  ck 


(3.32) 


□ 


Proposition  3.4  The  following  estimate  holds 


| \Q\ 1  |0,n 


<  C\\6o\\h 


where  C  depends  on  [  [ cro 1 1 s  for  s  >  {1  +n/2,/i}. 
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Remark.  With  a  few  necessary  changes,  the  proof  of  Proposition  3.3will  still  work.  Ob¬ 
serve  that  the  only  major  difference  to  the  analysis  before  is  that  each  one  of  the  transport 
equations  now  possesses  an  extra  term,  which  may  be  viewed  as  a  “source”  term.  Certainly, 
these  extra  terms  are  responsible  for  the  8a  term  in  the  estimate  of  our  main  theorem. 

We  may  also  rewrite  the  right  hand  side  of  (3.15)  in  terms  of  a^—ie9  plus  a  source  term 
containing  8a.  The  desired  estimate  then  follows  by  applying  Lemma  3.4,  together  with 
algebraic  estimates,  various  forms  of  Schauder’s  Lemma,  Rauch’s  Lemma  Proposition  2.1, 
and  commutator  lemmas. 

Note  that  the  transport  equations  of  a*  may  be  viewed  as  the  first  order  perturbation  of 
the  corresponding  equations  of  bk.  Hence  near  |x|  =  0,  by  the  main  assumption  (44)  8a  =  0, 
ak  —  0  for  k  =  1,  •  •  •,  /j. 

Next,  we  want  to  get  an  estimate  of  ||C2i^ig^it>|t=T||o,0'- 
Proposition  3.5  For  Qi,  w,  T .  and  0\  as  previously  defined,  we  have 

||Q,P,Tw|,=r||0,SV<  £11*11 

where  the  constant  C  depends  on  j | cr0 [  | s. ,  for  s  >  1  -f  n/2. 


Proof.  Since  iv  solves  the  dual  problem  (2.9),  which  has  smooth  right  hand  side,  we  need  a 
characteristic  trace  regularity  result  for  w.  Note  that  ES(Q i)  is  away  from  the  characteristic 
variety  of  the  wave  equation.  In  fact,  an  application  of  a  characteristic  trace  regularity  result, 
Corollary  2  in  [2],  yields  that  there  exists  a  <f>  €  C’^°(IRn+1),  such  that 

<  C'HHI,  <  C||*|| 

the  second  inequality  comes  from  the  standard  energy  estimates.  □ 

(b).  The  estimate  of  fij, 

Because  of  Proposition  3.4,  in  order  to  bound  i)^,  it  suffices  to  estimate  | \Q2fi1  ^^|t=r||o,n'- 
This  requires  the  microlocal  trace  regularity  result  Lemma  3.4  and  the  regularity  analysis  of 
the  dual  problem  Lemma  3.5.  The  estimate  can  be  easily  derived  by  applying  Lemma  3.4 
and  Lemma  3.5. 


Proposition  3.6  For  Q2,  w,  and  cq  defined  above.  The  following  holds 

On 

with  the  constant  C  depending  on  ||<7o||s  for  s  >  1  +??,/ 2 

(c).  The  estimate  of  /l  c 

We  want  to  estimate 

f  id 

/  Q-np-an  (f>\  —wds  . 

Ji=T{x)  p  On 
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Obviously,  this  term  has  to  be  handled  very  differently.  In  fact,  the  transport  operator  (the 
operator  appears  in  the  transport  equations)  is  elliptic  on  ES(Q i),  by  Lemma  3.1  the  L 2 
norm  of  Qifah  may  be  bounded  by  the  II1'  norm  of  Qifa0.  However,  there  is  no  ellipticity 
on  ES{Q 2),  and  the  only  direct  estimate  one  may  expect  is  by  applying  Lemma  3.2,  which 
implies  that  essentially  the  L2  norm  of  Q2dah  has  to  be  bounded  by  the  H21'  norm  of  </>a0. 

To  cure  this  difficulty,  we  introduce  a.  “bootstrap”  argument  by  taking  advantage  of  the 
transport  equations  and  the  fact  that  Q2d1~w\t=T  is  a.  much  smoother  term  (it  would  be 
smooth  if  (To  were  smooth).  Again,  we  start  with  the  corresponding  estimate  for  b q. 

Proposition  3.7  The  estimate 

I  [  Ch^-bh  (t>\~wds\  <  C||T|| 

Jt-T(x)  p  an 

holds  with  the  constant  C  depending  on  ||rr0||s>  for  s  >  l\  +  (n  —  3)/2. 


Proof.  Without  loss  of  generalities,  Q2  may  be  assumed  to  be  self-adjoint.  Thus 

d  .  f  1  ?  _  ,  d 

r{x)  p  Oil  Jt=r(x)  p 


f  \c)  f  1  d 

/  Q-20~hi,  Oj  -z—wds  =  /  -hi  <j)Q2<i>i^-ivds. 

Jt=r(x)  p  on  Jt=T(r)  p  on 


(3.33) 


Since  Q24>i-§^w  is  much  smoother  than  the  other  part,  one  wants  to  apply  transport  equations 
to  reduce  6/,  to  the  one  requires  the  least  regularity. 

Consider  a  problem 


(Vr-Vrilj  =  de^Q^—w)  |,=t(x) 


Hi 


0  near  lad  =  0  . 


(3.34) 

(3.35) 


Then  (3.33)  becomes 


1  ,  d 
d-lq 

H=T[X)  p 


,  4-wds  =  /  Vr  •  V6n  eqWi ds  . 

OH  Jt=r(3:) 


From  the  transport  equation  (3.16).  we  have 

/  VT-VbheqW,ds  -  [  (A(bi1^ieq)/2  —  |V^r|26/1_1e9/2  —  A</6/1_1e?/2 

Jt  —  T(x)  Jt=T(x) 

-V e")Vq0  +  bll_1eqVqWqo)W1ds 
—  I  b\  ieq  P\W\ds 

Jt  =  r{x) 

where  P  =  A  —  j  V <7 j 2  —  Aq  -f  V  •  Vcr0  -f  Vq  ■  Vq0. 

One  can  then  introduce  another  problem 

( Vr  •  V)*H'2  =  <j>PWi  (3.36) 

IT  2  =  0  near  |x|  =  0  .  (3.37) 
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Hence 


[  Q26-K  ©i  4-wds  =  /  Vr  •  S7b,^1eqPW2ds  . 

Jt-r(x)  p  On  Jt-r(x) 

We  can  continue  this  procedure.  Indeed,  it  follows,  after  [/i/2]  steps  where  [/i/2]  =  l\/2  if  l\ 
is  even  and  is  (l\  —  l)/2  if  l\  is  odd.  that 


f  =  f  Vt  ■  Vb[h/2]e9PW[i  /2]+ids 

Jt=T(x)  p  OH  Jt= T(X) 


where  W,  solves 


(Vr-VyWi  =  <f>PWi-i 

If7;  =  0  near  I  x  I  =  0 


(3.38) 

(3.39) 


for  i  =  1,  •  •  •,  [h/2]  +  1. 

Therefore,  by  the  Cauchy-Schwatz  inequality,  Proposition  2.1, 


and  Rauch’s  Lemma 


9  ,  I 

Wn 


< 


c  \\yb[hi2)tq  II/,  — 2[q/2]  |  \PW[h/2]+l  \\2[h/2]-h 


(3.40) 


with  C  depending  on  1 1 er0 1 1 ,  for  some  s  >  n/2.  Since  6,e9  solves  the  *-th  transport  equation, 
Lemma  3.2  then  implies  that 


\\tl>{ll/>]Ol\\ll-2{ll/2}  <  C' I  |0cro|  |/,  .  (3.41) 

Next,  applying  Lemma  3.2,  after  some  simple  calculations,  we  can  obtain 

rnw.lto./jw,  <  c\\^°-«QiM-2-w)\t=rM\h-i  ■ 

Now  Lemma  3.4  and  Lemma  3.5  may  be  applied  to  get 

\\JJW{h/2]+l\\2[l1/2]-h  <  C'PII 

where  the  constant  C  depends  on  1 1 <r0 1 1 s  with  s  >  lx  —  1/2  +  n/2.  □ 

Once  again,  a  variant  of  the  proof  of  Proposition  3.6,  with  no  additional  technical  diffi¬ 
culty,  will  lead  to  the  following  result. 

Proposition  3.8  The  estimate 

|  /  Qifh-o-i,  (pi—  ivds\  <  C'||^cr||/1||’I,|| 

Jt=r(x)  p  on 

holds  with  the  constant  C  depending  on  ||<t0||s,  for  s  >  l\  +  (n  —  l)/2. 
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